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Abstract 

We apply Renormalization Group techniques to the Real Time formulation of 
thermal field theory. Due to the separation between the T = and the T / parts 
of the propagator in this formalism, one can derive exact evolution equations for 
the Green functions describing the effect of integrating out thermal fluctuations of 
increasing wavelengths, the initial conditions being the renormalized Green functions 
of the T = theory. As a first application, we study the phase transition for the 
real scalar theory, computing the order of the transition, the critical temperature, 
and critical exponents, in different approximations to the evolution equations for the 
scalar potential. 
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1 Introduction 



The dynamics of a second order, or weakly first order, phase transition is governed by long 
wavelength fluctuations of the order parameter. These fluctuations, whose length scale is 
much larger than the inverse temperature of the system, are essentially classical, since the 
probability of a quantum fluctuation over such scales is highly suppressed. For this reason, 
the details of the microscopic theory are not relevant to the description of the critical 
behaviour of the system, which can be successfully studied by using classical models, such 
as the Ising model for ferromagnets or the Ginsburg-Landau theory for superconductors. 

In these models, a free energy can be defined as a functional of a macroscopic order 
parameter, which is allowed to vary only on large length scales. The effect of the short wave- 
length (quantum and thermal) fluctuations is incorporated by the parameters appearing in 
the free energy (masses, coupling constants, etc.), which are usually treated as phenomeno- 
logical parameters but should in principle be computable, starting from the underlying 
theory. 

The Wilson Renormalization Group (RG) |l| provides the natural framework in which 
this procedure can be systematically carried out. The main idea is to start from a micro- 
scopic theory, the parameters of which are supposed to be known, and then progressively 
integrate out the high frequency modes of the order parameter down to some infra-red 
cutoff A. In this way one obtains a coarse-grained order parameter and the corresponding 
effective action, which can be used as the relevant tool to describe the system. 

In this paper, we will apply this idea to quantum field theory at finite temperature. 
Our approach will be the following. First, we will assume that we "know" the zero- 
temperature renormalized quantum field theory, which means that there exists some reli- 
able approximation method (perturbation theory, lattice simulations, etc.) to compute the 
zero-temperature renormalized Green functions. In this way, the renormalization constants 
can be fixed by the experimental measurements and all the parameters of the theory are 
known. In a RG language, we assume that all the quantum fluctuations have already been 
integrated out. 

Second, we will integrate out thermal modes only, for frequencies higher than the in- 
frared cutoff A. For non-vanishing values for A, we will obtain coarse-grained order pa- 
rameter and free energy, which are the appropriate objects to study the dynamics of long 
wavelength thermal fluctuations. In the limit A — > 0, we will obtain the finite temperature 
quantum field theory in thermal equilibrium. 

In this approach, since all the quantum fluctuations are integrated out from the begin- 
ning, it will be possible to relate the thermal field theory and the physical (renormalized) 
quantum field theory at zero temperature in a transparent way. 

The ideal framework to perform this program of coarse-graining of thermal fluctuations 
is the Real Time (RT) formulation of thermal field theories ||, in which the thermal part in 
the free propagators is well separated from the zero-temperature quantistic one. As is well 
known, the price to pay for this is a doubling of the number of degrees of freedom. These 
"ghost" fields are necessary to cancel the so-called pinch singularities, which are due to the 



1 



fact that the thermal fluctuations are on-shell. Because of these technical complications 
Matsubara's Imaginary Time (IT) approach [[J is more popular in the literature. 

In particular, the coarse-graining procedure can be precisely formulated in the so-called 
Closed-Time-Path (CTP) formalism 0] (see refs. ||, [J for a clear presentation). Indeed, 
the description of a system in which only the short wavelength fluctuations are in thermal 
equilibrium can be achieved by modifying the density matrix with respect to the thermal 
one, and the CTP was designed just to describe systems with a generic density matrix. 

Anyway, the reader not familiar with the CTP formalism should not worry too much 
about it. As we will discuss in Appendix A, the modification of the density matrix that we 
will consider is equivalent to working in the more familiar RT formulation of Niemi and Se- 
menoff p], with a modified distribution function, given by the Bose-Einstein distribution 
function (or the Fermi-Dirac, for fermions) multiplied by the cutoff function. 

As a first stage of this program, in this paper we illustrate our method by considering 
the well-studied self-interacting real scalar theory. As is well known, this model belongs to 
the universality class of the Ising model, and then it has a second order phase transition. 
However, perturbation theory fails to reproduce this result even after the resummation of 
daisy and super-daisy diagrams ||, unless the gap equations are solved to 0(A 2 ), A being 
the quartic coupling constant ||. 

We derive the "exact" Wilson RG equation for the effective potential and approximate 
it by expanding in powers of space-time derivatives of the field. As we will discuss, even 
after our approximations, the effective potential computed in this approach includes more 
contributions than the super-daisy resummed one in H and the prediction of a second order 
phase transition will emerge in a clear way. In particular, we will see that the running of the 
coupling constant, which is neglected in perturbation theory, has a very strong effect. At 
the critical temperature, the effective coupling constant vanishes together with the thermal 
mass, and this cures the severe infra-red problems encountered in perturbation theory. 

The paper is organized as follows. In Section 2 we define the cutoff thermal field 
theory and discuss its relationship with other RG approaches. In Section 3 we derive 
the RG flow equations for the quantities of interest. In Section 4 the evolution equation 
for the cutoff effective potential is approximated in various ways in order to perform the 
numerical analysis. Section 5 contains the results of our numerical study and Section 6 
some conclusions and further comments. In Appendix A we briefly describe the Closed 
Time Path method and how it applies to our case. Appendix B contains a discussion of 
the problem of pinch singularities appearing in the kernel of RG flow equations. 



2 Cutoff thermal field theory 

As explained in the introduction, our aim is to define a coarse-graining procedure in which 
all the quantum fluctuations, and the high frequency thermal ones down to an infra-red 
cutoff A, are integrated out. In this section we will construct a path integral representation 
of the resulting generating functional of coarse-grained Green functions, in the case of a 
self-interacting real scalar field. 
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2.1 Cutoff propagators 



Let us consider a free field, expanded in the usual manner in terms of annihilation and 
creation operators 

4>{ x ) = J ^f2u k ^ exp (~ ik ' x ) + ^ ex P( ik ' x )] ' C 1 ) 

with ko = uj k = \Jk 2 + m 2 ; m 2 represents here the zero-temperature, physical mass of the 
scalar particle. 

In the case of an ideal gas in thermal equilibrium at a temperature T = 1/(3 one has || 
(a\a k :)p = {2-Kf2u k N{u k )8{k-k'), 

(2) 

(a k al) p = (2n) 3 2uj k [l + N(uj k )}5(k-k'), 

where N(ko) = [exp(/3£;o) — is the Bose-Einstein distribution function. The thermal 
averages in eqs. (0) are both temperature-dependent. In particular, in the second line, 
we have a "thermal" contribution, the N(uj k ) inside the square brackets, and a T = 0, 
"quantum" contribution, the "1". Of course, the commutator [afc,a|./] being a c-number, 
its expectation value is independent on the state and is the same as in the zero-temperature 
vacuum 

{[a k ,al,})p = (2Tif2u k 5(k-k') } (3) 

independent of T. We propose to modify the above relations by introducing an infrared 
cutoff on the thermal part only, that is 

(4a fc ,)£ = (27i) 3 2u k N(u k )Q(\k\,A)5(k-k'), 

(4) 

(a k al)% = (2n) 3 2u k il + N(u k )Q(\k\,A)]6(k-k>), 

where, in general, 0(|/c|, A) is a cutoff function which is 1 for \k\ > A and rapidly vanishes 
for \k\ < A. The simplest choice for 0(|fc|, A) is the Heavyside theta function 

9(|jfe|,A) -> 9(\k\ -A). (5) 

For simplicity, in the following we will always use the step function, even though in order to 
perform safe formal manipulations one should choose [[K|] a cutoff function which is always 



non- vanishing and of class C°°, and take only at the end the limit (|5|). 

The physical state described by the relations in (f|) is the one in which thermal equilib- 
rium is achieved only by the fast, high frequency, modes, from — ^ oo down to the cutoff 
scale \k\ ~ A. The corresponding density matrix is given by 



p = cexp 



d 3 k f3 kt A %a k 



(6) 



where (3 k ^ = f3uj k for \k\ > A and goes to infinity for |A;| < A. The quantum relation 

(K4]>£= (2n) 3 2u k 5(k-k'), (7) 
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analogous to eq. (H), is valid, in this state, for any value of the momenta, irrespective of 
the cutoff. Note that this would have not been the case if we had multiplied the whole 
right-hand side of the second eq. (pj) for the cutoff function G(|/c|, A). 

Now the cutoff Green functions are defined in the usual way, as statistical averages of 
Heisenberg fields ordered along a path C in the complex time plane |5], ||, running from 
— oo to +00 along the real time axis and back from +00 — ie to —00 — ie infinitesimally 
below it. As we will show in Appendix A, the effect of the non-thermal density matrix 
(^) can be accounted for by simply working in the usual equilibrium formalism, but with 
the Bose-Einstein distribution function multiplied by the cutoff function. For instance, the 
two-point Green function is 

G®(x, x') = 9 c (t - tf)C>{x, x') + 6 c (t' - t)C<{x, x'), (8) 

where 9 c (t—t') is the generalization of the theta function on the contour C and the two-point 
correlation functions are 

C>{x, x') = ($(x)$(x'))% = C<(x\ x) . (9) 

The free cutoff propagator D^\x — x') = — iG^\x, x') can be computed in the standard 
way, by substituting (D) in (||) and (||), and using (|). One obtains 



d 4 k 

x [6(ko)9 c (t - t') + 6(-k )6 c (t f -t) + N(\k \,A)] , (10) 



iD { l\x - x') = J tS^A; 2 - m 2 )e- ik < x -^ 



where 

N(\k \,A) = N(\k \) 8(\k\ - A) . (11) 
The propagator can be seen as a 2 x 2 matrix whose components are given by 



(12) 



where t, t' lie on the real time axis and we have omitted the spatial arguments of the 
propagator. 





-t') = 


Df{t- 








~ t') = 




— ie) 


-(f-ie)), 




-*') = 


Df{t- 


-(*- 


ie)), 


Df\t 


~t') = 




— ie) 


-t'), 



Using expression (10) and then Fourier transforming, we obtain the RT cutoff propa- 
gator in momentum space 



D A (k) 
where 



/ A (A - A*)9(-k ) 

{ (A - A*)6(k ) -A* 



+ {A -A*)N{\k \,A)B, (13) 



A o = T2 W- ( 14 ) 

k z — m z + ie 
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and 



/ 1 1 



B= . (15) 

Note the separation between the T = and the thermal part. In the e — > limit, we have 

A -A* — ► -2in5(k 2 -m 2 ), (16) 

that is, the thermal part of the tree level propagator involves on-shell degrees of freedom 
only (only real particles belong to the thermal bath). 

Since it will be useful in the following, we give here also the expression for the derivative 
of the propagator with respect to the cutoff 



^D A (k) = 2m6(k 2 - m 2 ) 5{\k\ - A) N(\k \) B 



[IT) 



and for the inverse of the propagator 

1 / 

D A {ky l 



A* (A - A*)6(-k ) 



A AS 
A -A 



a as 



V (A - A*)9(k 
S JV(|fco|,A) 



-An 




(18) 



2.2 Cutoff effective action 

By introducing independent sources ji and j% for the two pieces of the contour C (the real 
time axis and the t — ie one, respectively), the path integral representation of the generating 
functional of RT cutoff Green functions can be written as @ 

Z A [j] = J [#!][# 2 ] exp^ {±tr • Dl 1 • + S[<f>] + tr j • <f>] , (19) 

where the trace means integration over momenta while the dot represents the sum over 
field indices 



a,b=X,2 ' 



(20) 



trj-0 = ^ J J^uJa(-p)Mp) 



a=l,2 



(27T)' 

and >S[0] is the bare interaction action 

S[0] = S[0 1 ]-S[0 2 ]. (21) 

Notice that in eq. ([T^) it is assumed that the usual procedure of renormalization of the 
perturbative ultraviolet divergences of the zero-temperature theory has been carried out. 
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Namely we assumed a regulator, whose precise form is irrelevant for our discussion, and a 
set of zero-temperature renormalization conditions. 



In the usual interpretation of the RT formalism, X and 02 are respectively the "physical" 
field and the "ghost" field. 

From the expression of the propagator in ( O ) one can immediately realize that, in the 
limit A — > oo (which in practice, due to the exponential behaviour of the Bose-Einstein 
distribution function, means A ^> T), none of the thermal modes are propagating inside the 
loops. In this case by taking the functional derivatives of (|19|) with respect to the source 
ji one obtains exactly the zero-temperature, fully renormalized, quantum field theory. On 
the other hand, in the opposite limit A — > the propagator tends to the equilibrium, Real 
Time, finite temperature propagator, and consequently the generating functional in (|19|) 
gives the full finite temperature theory in thermal equilibrium. 

We define as usual the cutoff effective action as the Legendre transform of the generating 
functional of the connected Green functions W\[j] = — i log Z A [j] |L1|, [12], p~3|| : 



lti<p-Di 1 -<p + r A [<f ) } = w A [j]-tij 



(22) 



where we have isolated the free part of the cutoff effective action and used for the classical 
fields the same notation as for the quantum fields. 

As discussed in [|], the free energy of the system is given by the functional f A [<£>], defined 
by the relation 

5T A [<p] 5I\[0] 



5ip 



(23) 



!>i=<y>2=¥> 



The tadpole F A \ip), which will play an important role in the following, is found by evalu- 
ating (|23| ) for a space-time independent field configuration 



f(i) 

1 A 



STaM 



5ip 



(24) 



i(3=const. 



We conclude this subsection by recalling an important property. By inspection, one finds 
that the cutoff effective action has a discrete Z2 symmetry p| 



r A [0i,0 2 ] = -r; 

This relation has the following consequences on the derivatives of the tadpole 



dip 



P2=<P 



dip 2 



E 

a,b=l,2 



5 3 r A [0] 



(25) 

(26) 
(27) 



where again ip is space-time independent. We will use these relations below. 
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2.3 Relation to the other formulations of the Wilson RG 



Before proceeding with the discussion, we would like to comment on the relationship be- 
tween our approach and the other formulations of the Wilson RG existing in the literature. 

The first application of the continuum Wilson RG philosophy to quantum field theory at 
T = is due to Polchinski [|n| . The motivation of that work was to provide a demonstration 



of perturbative renormalizability without resorting to diagrammatic techniques. 

In this formulation, which we will indicate in the following as the Polchinski RG, only the 
modes between an infra-red cutoff A and an ultraviolet cutoff Ao are allowed to propagate 
inside loops. The bare Euclidean Lagrangian depends only on the UV scale A , which 
is eventually sent to infinity, and the renormalization group flow describes the effect of 
quantum fluctuations coming into play as the infra-red cutoff A is lowered. This can be 
achieved in a way formally analogous to what we have done in this section, the main 
difference being in the form of the cutoff propagator. Namely, in this case the cutoff 
procedure has to be applied to the full propagator, and not only to a part of it, as we have 
done here. That is 

D AAo (k) = e(k, A)6(A , k)D(k) , (28) 

where D(k) is the tree level propagator. The internal legs of the Green functions generated 
by the analogues of our Z\, W\, and T\, now carry only momenta A < k < A . Therefore 
the 1PI Green functions are given by the bare couplings for A = A (since no modes are 
propagating inside loops) and by the fully renormalized ones for A = (since in this limit 
the quantum corrections have been included at any momentum scale). 

On the other hand, since the cutoff procedure described in this section modifies the 
thermal part of the propagator only, the Green functions defined in this paper contain the 
full quantum corrections, independently of the value of A. Lowering A, we are introducing 
new thermal modes only. 

In short, the Polchinski formulation of the Wilson RG interpolates between the bare 
and the physical (renormalized) theory, whereas the formulation presented in this paper 
interpolates between the physical theory at T = and the physical theory at T ^ 0. 

Since the two RG's describe two different physical problems, also the boundary con- 
ditions will be different. As stated previously, at A = Ao in the Polchinski RG we have 
the bare Lagrangian. Thus, one should impose two different classes of boundary conditions 
|T(]| , |T2|| . The first one involves all the "irrelevant" operators, the ones with dimension larger 
than 4. These should satisfy the power counting in A at the ultraviolet scale A = A . The 
remaining, "relevant" , couplings have to be fixed by the matching of the renormalized the- 
ory and the physical observables, and so the initial conditions for them should be given at 
A = 0. 

On the other hand, in the approach presented in this paper, the physical theory corre- 
sponds to the initial condition at A — ► oo, so that all the couplings, both the relevant and 
the irrelevant ones, have to be fixed at this point. 

The Wilson RG method has already been applied to the study of the T ^ quan- 
tum field theory in the IT formalism in ref. [16| and more recently in [17|. The thermal 
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propagator in the IT formalism [18] is obtained from the Euclidean propagator after the 

(29) 



replacement k —> u n , where u n = 2nTn are the Matsubara frequencies 

1 1 



k 2 + m 2 (2vrTn) 2 + \k\ 2 + m 2 

and the ko integration is consistently replaced by TJ2 n - 

While in the RT formalism the quantum and the thermal part of the propagator can be 
clearly identified, as in fll3"D , in the Imaginary Time formalism this separation cannot be 
achieved, as we read from (p9|). So, in order to formulate a Wilson RG in this case, we can 
only modify the thermal propagator as in (p8|) , imposing an IR cutoff on the combination 
(2nTn) 2 + \k\ 2 . 

In this case we recover the bare theory for A = A , since no modes, either quantum or 
thermal, have been integrated in the loops, while for A = we have the physical theory at 
T ^ 0. This formulation of the RG describes the effect of thermal and quantum fluctua- 
tions at the same time, and it is different from both Polchinski's and ours. Concerning the 
boundary conditions in this case, it should be noticed that since the RG now interpolates 
between the bare theory and the T / physical theory, there is no value of A that cor- 
responds to the physical theory at T = 0, the one that is supposed to be directly related 
to the measurable observables. Then, in order to match the running parameters with the 
physical theory, a preliminary step has to be performed, consisting in the derivation of the 
bare parameters from the renormalized, T = theory |TB|, |T7|. 



3 Renormalization group flow equations 



In this section we study the A-dependence of the cutoff effective action. By taking the 
derivative with respect to A of (|I9| ) we obtain the evolution equation for Z A [j] 



A^-Z A \j] = -- tr — ■ A^-D^ 1 ■ 5 -ZM . 
dA m 2 5j dA A Sj Ui 



(30) 



Choosing as initial conditions for Z A [j] at A ^> T the full renormalized theory at zero tem- 
perature, the above evolution equation describes the effect of the inclusion of the thermal 
fluctuations at the momentum scale \k\ = A. 

The evolution equation for the generating functional of the connected Green functions 
W\[j] = — i log Z\[j] can be derived straightforwardly from eq. ( |30"D 



-tr 



\d_ n -i 
dA ' SjSj 



2 Si OA A 



(31) 



By using this equation in 
found: 



Sj OA Sj 
the flow equation for the cutoff effective action T^[4>] is 



i 

2 tr 



D A l + 



s 2 r A \ 



(32) 
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Notice that the terms coming from the A-dependence of the classical field <fi cancel out, 
which is a well-known general property of the Legendre transform. The evolution equations 
for the various vertices can be found by expanding the r.h.s. in powers of <fi |12, |i~3|j . 

Deriving eq. (|32|) with respect to 4>i an d then setting <pi and 2 equal to a constant 
background (p, we obtain the evolution equation for the tadpole, defined in eq. (0). 



^A 1 + Sa(^) 



-1 



D^ + Sa^) *-rf{<p)\, (33) 



5r A [0] 



where 

K OCpi OCpj 0<p X 01=02=¥ , 

and we introduced the field-dependent self-energy matrix 



(2v)W\p + p>) [E A (p; 



(34) 



(35) 



'1=<P2=V 



First of all, we need to study the kernel of the evolution equation for the tadpole, that is 

-l . d _ 1 r_ _ , .i-l 



K A (k; if) = -% Dl v + £ A (¥>) 



(36) 



^a 1 + M<P) 



-i 



. This can be 



The kernel contains the "full" cutoff matrix propagator 
obtained (by assuming a Schwinger-Dyson equation 0) from eq. (|13|) by substituting A 
with 

A^A g o (-n t A )- ^ 2 _ m2 + nA(fe;y)+fe ,37) 

where the quantity n A (A;; <£>) is related to the 11 component of the self-energy matrix by 

ReIT A (£;; (p) = Re [E A (fc; <p)] n 

1 (38) 



Imn A (fc; V 5 ) 



Im [E A (fc; V 5 )] 



ii 



l + 2JV(|feo|,A) 

Collecting the above formulae we can compute the kernel, which turns out to be 



K A (k; <p) = -i(Ao - A* 



A<J(|fc|-A)iV(|fco|)fl 



(39) 



Notice that if one uses a cutoff function different from the step function, the kernel is 
obtained from (|39"D by substituting S(\k\ — A) with the A-derivative of the new cutoff 
function. 

Comparing (^) with the derivative with respect to A of the tree level propagator, which 
we have computed in fllTD, we see that the delta function which forces the momenta on the 
zero temperature mass-shell has been replaced by 



-2np A (k; <p) + *A A*(II A - IT A )(1 + R A + R* A + R A R 



A7> 



(40) 



where 



p A (k; V ) = ~[A A -A\] (41) 



is the interacting spectral function and 

oo 



n=l 



The second contribution to the kernel in eq. (flCf) exhibits pinch singularities of the form 
A™ Aq™. It is also proportional to the imaginary part of the self-energy ILa(A;), which 
is non-vanishing along the whole real axis in the k complex planeFl. As will be shown in 

(3) 

Appendix B, these pinch singularities cancel with analogous contributions coming from T A 



once the kernel ( f40|) is inserted in the evolution equation for the tadpole to give 

a m^ 1)( ^ = "4 / (0 ^ - a)(a ° - a ° } ^a| N{lk ° l)T ® {k] ~ k] ^ ' (42) 



where 



(here the trace is over 1, 2 field indices) and 1^ nas been defined in (0). 

It is important to stress here that the above equation is exact, since no approximation 
such as perturbative expansion or truncation has been performed up to now. 

In the rest of this section we will discuss the evolution equation obtained neglecting 
the imaginary part of the self energy on-shell, which in perturbation theory arises only at 
two- loops [13]. We will use this approximation as the starting point for the numerical 



analysis which we will describe in sect. 4. In this approximation, the kernel in ( [40]) can be 
written as 

\ \ (44) 

k 2 - m 2 + ReU A (k; <p)+ie k 2 - m 2 + R e n A (fc; tp) - is 

so that the only singularities are the poles in the ko complex plane satisfying k 2 , = \k\ 2 + 
m 2 — ReIlA(fc; (p) ±ie = 0, and it vanishes in the rest of the complex plane. Notice that 
the pinch singularities have disappeared in this approximation. 

In the case in which \k\ 2 + m 2 — KeU A (k; </?) > the poles lie infinitesimally close 
to the real axis, and ([44|) reduces to 2-nbik 2 — m 2 + KeU A (k; ip)). On the other hand, 
if \k\ 2 + m 2 — ReIlA(A;; p) < 0, as can be the case when there is spontaneous symmetry 
breaking, the poles lie on the imaginary axis, and the contributions of the two pieces in 
(P|) to the integration along the real k axis cancel. 



Then, the quantity ([40D appearing in the kernel reduces, in the e — > limit, to 

2(A - AS) ^| - 2n5(k 2 -m 2 + ReU A (k; <p)) 9(\k\ 2 + m 2 - ReU A (k; <p)) . (45) 

The physical meaning of this kernel is straightforward. As the infrared cutoff A is lowered, 
the mass associated to the new modes coming into thermal equilibrium is not given by 



1 We thank S. Jeon for drawing our attention on this point. 
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the T = pole mass m 2 , but by the thermal mass, that is the pole of the full propagator 
obtained by integrating over the high momentum modes (\k\ > A), already in thermal 
equilibrium. This comes out quite naturally in this formalism, whereas in perturbation 
theory one has to perform ad hoc resummations in order to cure the infrared divergences. 
We will come back to the comparison between this approach and perturbation theory in 
the following. 

The other interesting feature of this kernel has to do with the case of spontaneous 
symmetry breaking. In perturbation theory, one is not allowed to compute the effective 
potential for values of the background field close to the symmetric phase, if the temperature 
is less than the critical temperature of the phase transition. This is because the thermal 
mass squared is negative, in this region, which gives rise to a complex effective potential. 
As showed by Weinberg and Wu |I9[ , the real part of this effective potential can still be 
interpreted as the energy density of a spatially homogeneous, although unstable, state. 

On the other hand, in this approach, we see that the only modes giving rise to a thermal 
evolution are those with a real energy, and the running stops as soon as the energy squared 
becomes negative, so that no imaginary parts for the effective potential are generated. 



Using (f45p , the evolution equation for the tadpole (|^) now takes the simple form 



(2nY 

x6(\k\ 2 + m 2 -ReU A (k; (p)) (jfe; -Jfe; <p) , (46) 

which will be studied numerically in the next section in an approximation scheme based 
on a derivative expansion and truncations. 

Before concluding this section, we would like to briefly discuss the infra-red limit of the 
evolution equation for the tadpole ([46]). As the energy of the thermal modes becomes much 
smaller than T, the Bose-Einstein distribution function can be approximated as 



N(\k \)^-- + . (47) 




Integrating in k in (^) and keeping only the first term in the expansion for N(\ko\) in 
(fl7j) the evolution equation for the tadpole becomes 

9 pw/.^, AT f d " k ™a Tf(k--k-v) 



if 



Equation (fig) shows that, in the limit of eq. fl49|), and neglecting the imaginary part of the 
self-energy, the exact four-dimensional running at finite temperature reduces to a purely 
three-dimensional one, at zero temperature. The momentum integral is now a genuine 
three-dimensional loop integral with the zero-temperature propagator. Equation (48) is, 



a— rr (vo * ~=- / 7^ m - a) r A ; p ; r ; r - (4* 

c*A 2 J (2-k) 6 \k\ 2 + m 2 + ReIl A (k; (p) 



lu 2 a = A 2 + Ren A (fc = w A> |^| = A; <p) < T 2 . (49) 



in fact, the same equation as would have been obtained in the Polchinski RG in three 
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dimensions and T = 0, provided the following matching between the three-dimensional 
and the four-dimensional couplings had been done: 

ri 1} (^D=3 = 

ReIl A (k; (f) D =3 = ReU A (k = u k ,k; tp) D=A (50) 

Ia (£; —k; <p) D =3 = Vrf a } (^o = Uk, fc; &o = ^k, -fc; v)d=4 • 

Obviously, the above matching has to be performed at a value of A such that the condition 
(HH) is fulfilled, that is in the three-dimensional regime. So, the question arises whether we 
can neglect, in a first approximation, the four- dimensional running from A = Ao 3> T to 
some value of A, A 3 £>, which is inside the three-dimensional regime. The answer is clearly 
negative, as we can read from Fig. [TJ. 
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Figure 1: The running of the thermal mass (Vj[((p = 0), here given in units of v 2 , v being the 
zero temperature vev) with respect to the ratio A/T. The region to the left of the dashed 
line corresponds to the three-dimensional regime (see text). We have fixed T = 1.2 T c and 
A = 0.1. 

In this picture, we have plotted the value of the thermal mass in (p = (see Sect. 5) as 
a function of the ratio A/T. The area to the left of the dashed line represents the region 
in which the Bose-Einstein distribution function can be approximated by T/\k Q \ with an 
accuracy better than 10%, and consequently the three-dimensional RG equation (fP|) is 
valid to that accuracy. 

As we can see, most of the running of the mass, from the zero-temperature value 
(corresponding to A/T ^> 1) to the finite temperature one, takes place before the three- 
dimensional region is reached. In other words, in order to match the three-dimensional 
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couplings with the physical four- dimensional ones, it is essential to keep the pure four- 
dimensional running (at under control, and the study of the pure three-dimensional 
running is, in general, not justified. 

On the other hand, dimensional reduction is justified when one is interested in the study 
of the critical theory. By universality, we know that this theory is insensitive to the "initial 
conditions" of the RG equations, so that the matching with the physical theory at zero 
temperature is not very important for the computation of universal quantities. 

In our numerical analysis we will always consider the four- dimensional RG equations 
derived in the previous sections, without dimensional reduction. This approach will allow 
us to study both the universal quantities (critical exponents) and the non-universal ones 
(critical temperature, thermal masses), and to keep the transition from the four- dimensional 
to the effectively three-dimensional critical theory under control. At the end of sect. 5 we 
will comment on the possibility of using a dimensionally-reduced form of eq. fl42|) in order 
to study critical regime of the theory. 



Derivative expansion, truncations, and comparison 
with perturbation theory 



An analytic solution of the exact evolution equation for the tadpole in eq. fl42|) is not avail- 
able. In principle one has to know the momentum dependence of both the full self-energy 
and the vertex f , i.e. of the first and second derivative of the tadpole with respect to the 
field (p. A similar problem is encountered in the applications of the Polchinski RG in the 
zero-temperature field theory. An approximation scheme based on a derivative expansion 
(or momentum-scale expansion pOJ) has proved to be very efficient in that context, at least 
in the case of the scalar theory. 

In order to perform a systematic derivative expansion in this case, we expand the free 
energy functional defined in ( p3|) in derivatives of the field <p(x) as follows: 

d A x 



imV - V A {<p) + \{d V ) 2 Z K { V ) + \{u ■ d V ) 2 Y K ( V ) + ■■■] (51) 



where the dots indicate higher derivative terms. Note the term containing the quadrivector 
of the thermal bath, u^, which has to be introduced for relativistic covariance. 



Stopping the expansion in ( pl[ ) at some order, plugging it into eq. (f±2|), and equating 
the coefficients of the terms with the same powers of derivatives of the field, we obtain 
evolution equations for the functions V\, Z^, Y\, and higher orders. 



A remark is in order. Usually in the Polchinski RG the expansion ( pT|) is not well 
defined if one uses a sharp cutoff function. The problem arises because in the zero external 
momentum limit the analogue of our eq. ( j3^) and its momentum derivatives contain ill- 
defined products of S(\k\ — A), coming from the derivative of the cut-off propagator, and 
0(\k\ — A), coming from the underived propagator [20|. Fortunately, this is not the case in 



our approach. The point is, again, that the theta function here appears only in the thermal 
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part of the propagator, and all the above mentioned products of theta and delta functions 
cancel one another by the same mechanism that ensures the cancellation of the "pinch" 
singularities in the RT formalism. As a consequence, the sharp cutoff limit can safely be 
taken even at higher orders in the derivative expansion. 



At the lowest order, the derivative expansion corresponds to neglecting the momentum 
dependence of the self-energy and the three-point vertex appearing in the r.h.s. of (f42|) . 
In this approximation we can use (|26|) in (|37]) and Aa can be written as Aa = [k 2 — 
V^tp)]' 1 , where the prime over the effective potential indicates derivation with respect 
to the field. Analogously the three-point vertex in ([43]) is equal to — V^"(y?). Note that 
since the imaginary part of the self-energy vanishes for zero external momenta, it does not 
contribute to the zero order in the derivative expansion. We get 



A k v ^ 



A 3 N(u A (cp)) 



47T 2 UJ A ((f) 



V^{p)9{A 2 + V'l{p)) 



(52) 



where 



Using the equilibrium Bose-Einstein distribution for N and integrating in (p we find the 
following evolution equation for the effective potential 



4a^ 



^ A 3 , 
- T 2^ l0g 



exp (-Pyftf + Vj^tp) 



9(A 2 + V^p)). 



(53) 



The same equation was found in the Matsubara formalism in [pT[ . 



Looking at the effective potential as a function of the two variables A and ip we see 
that the evolution equation is a non-linear partial derivative differential equation, with the 
initial condition V A=oa (p) being the renormalized effective potential at zero temperature. 



In principle, one could seek for a numerical solution of the evolution equation 
An alternative procedure, which will be followed in the present paper, is to take further 
derivatives of eq. (|52|) with respect to the field (p, so that the partial derivative differential 
equation is turned into a infinite system of ordinary first order differential equations in A 

, with 99-dependent coefficients, 



for the unknowns V A , V A , 



V" 



A k v ^ 



A k v ^ 



A 3 N{u A ) 



Air 2 



u A 



VM9(A 2 + VM) 



A 3 (N(u A ) 



'Ait 2 \ 



u A 



N(co A ) 
u A 



Km 2 e(A 2 + vM) 



-a^a mm] 



(54) 

In the following we will solve this system of equations at different orders of truncations, in 
order to test the reliability of this further approximation (besides the derivative expansion). 
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This also allows us to make a comparison between the present approach and the per 
turbative one. The 1-loop perturbative result and the daisy and super-daisy P, p|, p3 



resummed ones correspond to different approximations of the above system, truncated at 
most to second order, i.e. without taking into account the evolution of the third and higher 
order couplings. More precisely: 

(i) 1-loop perturbation theory corresponds to a truncation to the first equation 
only, the one for the tadpole, in which Vj[ and V^' are approximated by their tree level 
values: 

V^)-™%), TO^V; (55) 

(ii) Daisy resummation improves the 1-loop result by partially taking into account 
the evolution of the "mass" Vj[, replacing m 2 (ip) with a A-independent "thermal mass" 
in the evolution equation for the tadpole. The trilinear and quadrilinear couplings are 
approximated by the tree level ones, that is 



oo 







(56) 



where the function defined in (|54|), is computed by approximating V^(cp) with the tree 
level mass. 

(iii) Super-daisy resummation corresponds to an improvement of the previous ap- 
proximation consisting in using as the value of the thermal mass the solution of the 
"gap equation" [] 

Vj&<p) ~ ff%(<p) = m 2 (tp) + / dA'^A' rh 2 T {(p) 

Jo (57) 

VR{tp) a Xip VT(<P) * A. 

As we see, the (resummed) perturbative results do not take into account the evolution 
equations for the trilinear and higher orders couplings. Moreover, the A dependence of the 
second derivative of the effective potential is not considered, while a A-independent thermal 
mass is introduced. 

In the next section, we will see how these effects are in general very important, mainly 
at temperatures close to the critical one. In particular, the inclusion of the running up to 
the fourth derivative of the potential, will change the transition from a (weakly) first-order 
one to a second-order one. 



5 Numerical results 

We have solved the system in eq. (|5^ ) at different orders of truncation and for different 
values of (p. In this way we are able to reconstruct the shape of the tadpole, and then of the 

2 Different resummation schemes have been presented in the literature, which however agree with one 
another up to 0(/3), where (3 = XT/rriT is the effective expansion parameter. 
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effective potential. Since we are mainly interested in the study of the phase transition, we 
will solve the evolution equations for <p close to the origin, and running from A = A ^> T 
(in practice A > 10T will be enough, see Fig. [I]) down to A = 0. As initial condition in 
A = A we use the 1-loop effective potential at T = 0, 



rr ( s 1 2 2 . A 4 . ^(l) 



log 



m 2 ((p) 
-2/i 2 



(5? 



where 



2 / \ 2 ^ 

m {(f) = —fi + ~(p 



and the potential has been renormalized in the MS scheme with Q 2 = —2/i 2 . In principle, 
different approximation schemes to the T = effective potential could also be used, as given 
for instance by the Polchinski RG or by lattice computations. However, due to universality, 
the results for the critical theory should be almost insensitive to the details of the initial 
conditions, and the 1-loop approximation will be enough for our present purposes. 




Figure 2: The tadpole VJ( =0 (in units of v 3 ) as a function of (p, obtained by truncating 
to order II. The temperature has been chosen such that the effective potential has two 
degenerate minima. We have fixed v = 246 GeV and A = 0.1. 

In Fig. H] we show the tadpole V^ =0 as a function of cp at the critical temperature, 
when the system has been truncated to the second line. As we see the potential has three 
stationary points (the zeros of V^ =0 ), i.e. there are two minima, so that it describes a 
(weakly) first-order phase transition. In this case, we have fixed the temperature in such 
a way that the two minima are degenerate. This result agrees with the one obtained in 
1-loop resummed perturbation theory, see for instance ref. ||. 

However as soon as we turn the evolution for the third and fourth coupling on, as 
we did in Fig. |3|, we see that things change dramatically. In this case the tadpole has 
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Figure 3: The tadpole V^ =0 (in units of v 3 ) as a function of <p, obtained by truncating to 
order IV (continuous line) and VI (dash-dotted line). The temperature corresponds to the 
vanishing of the second derivative at the origin. We have fixed v = 246 GeV and A = 0.1. 

only one zero in (p — for temperatures higher than a critical value and two zeros for 
lower temperatures. There is no temperature interval in which a third zero is present, 
in other words the effective potential develops no barrier between the symmetric and the 
asymmetric phases. The critical temperature T c has been determined by requiring that the 
second derivative of the potential vanished. So, as soon as the running of the coupling is 
taken into account, the transition turns out to be of second-order, as would be expected 
from universality. The continuous and dash-dotted lines in this plot represent the results 
obtained by truncating to the fourth and sixth orders, respectively. 

The running of the coupling is indeed a strong effect, in particular close to the critical 
temperature, as we read from Fig. |j. The T = value for the coupling in this picture is 
chosen to be 0.1, and we see that, while for T far from T c the running is an ~ 10% effect, 
for T — > T c the coupling vanishes. 

The expansion parameter of the super-daisy resummation is given by (3 = XT/rriT, 
where my is the mass scale of the theory, in this case the thermal mass V^ =0 . As T — > T c 
the thermal mass vanishes, and (3 diverges. This is the source of the wild infra-red problems 
of resummed perturbation theory, which prevent computations for temperatures close to 
the critical one. Including the running of the coupling, which is also vanishing at T c , (3 
tends to a finite value, and the infrared problems are, if not cured, at least domesticated. 
This observation was already made in ref. [|I~6|j . 

Looking at the evolution equations, we see that the range of A for which there is a 
sizeable running of the parameters is determined by the function N(u\)/uj\. For A ^> T 
the distribution function is exponentially suppressed and the running is negligible, as we 
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Figure 4: The running of the coupling constant V£l Q ((p = 0) as a function of the tempera- 
ture, for T > T c . We have fixed v = 246 GeV and A = 0.1. 



read from Fig. [|. High energy modes are Boltzmann suppressed. 

On the other hand, for small A, the function goes approximately as T/(A 2 + V^'), and 
so there is running as long as A 2 > V£. We can then define a correlation length as the 
inverse of A c , such that 

Ac = V'L ■ (59) 
Since there are no thermal modes at wavelengths larger than the correlation length, inte- 
grating over A < A c gives a negligible effect. 

At the critical temperature we have V^' =0 = 0, and the correlation length diverges. We 
can define critical indices in the usual way (see for instance ref. |2^|) and see how their 
computation is improved with respect to perturbation theory in the present approach. 

At the lowest order in the derivative expansion, we have computed the critical exponent 
governing the scaling of the renormalized mass near the critical temperature, 

m A=0 ~ |T - T c \" , (60) 

and S, describing how the magnetization at the critical temperature, M(h, T = T c ) = 0, 
scales with the magnetic field h — V^, 

M(h, T = T c ) = ~ l/il 1 /* = (ytf/ & . (61) 



In Table [l] we summarize our results for the order of the phase transition, the critical 
temperature, and the critical exponents v and 5 obtained at the lowest order in the deriva- 
tive expansion, and truncating the system in eq. fl54|) at second, fourth, and sixth order 
respectively 
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Truncation 


A 


Tjv 


Order trans. 


V 


5 


V 


II 


0.1 


1.99 


I 








IV 


0.1 


2.01 


II 


0.53 


3.27 




IV+Z„ 


0.1 


2.02 


II 


0.53 


3.27 


0.015 


VI 


0.1 


2.02 


II 


0.58 


3.57 




Best values |24|, |25|| 






II 


0.63 


4.82 


0.032 



Table 1: Results for the order of the phase transition, the critical temperature, and the crit- 
ical exponents u, 5, and 77 at different orders of approximations to the evolution equations 
(see text). The last row contains the best values in the literature. 



In the last row we have listed the best results for the critical exponents obtained in the 
literature |23, |25| . We see that passing from the fourth to the sixth order of truncation we 



have a sizeable improvement in the result for v and S. 

Including the corrections of O ((dp) 2 ), keeping also the third and fourth terms in the 
expansion in (|5T|) , we can compute the critical exponent rj, which describes the scaling of 
the two-point renormalized Green function with respect to A for A — > at T = T c : 

ri 2) (0)=mi(l + Z A )~A 2 -" . (62) 

Taking again derivatives with respect to the field tp, we obtain a system of first order 
differential equations for Vf[, ■ ■ •, Za, Z' a , . . ., Ya, Y a , . . .. We have truncated this 
system at the fourth derivative for V\ and at the second one for Z/^ and Y\. Due to the 
Z2 symmetry of the model, and using the fact that the initial condition for Ya (which 
corresponds to its zero-temperature value) is zero, the effect of Ya and its derivatives is 
negligible for values of <p close to the origin, and can be neglected with respect to that of 
Z?y. Moreover, we have still neglected the imaginary part of the self-energy on-shell, which 
is only a two-loop effect in standard perturbation theory ]TJ|, 



The result is reported in the fourth row of Table [I]. As we can read, the effect of wave 
function renormalization is very small, as may be expected since in perturbation theory 
it arises at two loops. So, in this model, the derivative expansion turns out to be a very 
efficient approximation scheme. 

The accuracy of our results for the critical exponents is generally a lot better than 
the one obtained in other approaches to the four-dimensional theory, such as perturbation 



theory or the IT formulation of the RG [16|, both predicting v = 0.5. 



Nevertheless, the best way to compute the critical exponents remains the study of the 
three-dimensional critical theory (see refs. |24], ^] and, for a RG computation, [^(J ). Once 



one has established that the effective potential solution of eq. (|53| ) describes a second order 
phase transition, the universal quantities (the critical exponents) can be obtained studying 
the evolution around the fixed points of the critical (T = T c ) theory at large length-scales 
(A — > 0). For A < T c one expects that the theory becomes scale invariant. This is what 
indeed happens, as can be seen by rescaling V\ and tp according to Vx = A 3 T C V\ and 
p = VATc (p (notice that the temperature enters as expected by dimensional reduction 
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arguments). In terms of the new variables, eq. (|53D becomes 



3V A 



2n- 



log 



1 — exp ( — e 



t = log 



A 



As t — ► oo the explicit t-dependence in the RHS disappears, and we obtain the scale- 
invariant equation (dropping tildes and changing variables to absorb the factor in front of 
the logarithm) 

JV A + yv A - 3v A = log (i + v£) . (63) 

This is the RG flow equation for the zero temperature potential in 3 space-time dimensions 
at the lowest order in the derivative expansion and for a sharp momentum cutoff [27]. We 
see than that the critical behaviour of ( |53| ) is effectively described by a three dimensional 
T = theory. The fixed points and the corresponding critical indices of eq. ( |63"D can be 
computed numerically without truncations and the result for v is (see [B7| for details) 



v = 0.6895. 



On the other hand, if one is interested in studying the theory out of the critical regime, 
or in the case of a first order phase transition, the three-dimensional approach is no longer 
suited, and the full four-dimensional theory has to be addressed. In this case, as we have 
discussed at the end of sect. 3, the present formulation of the Real Time RG provides a 
clear connection between the theory at finite temperature and the theory at T = 0, which 
we are supposed to test in the laboratory. 



6 Conclusions and outlook 



In thermal field theory several difficulties are encountered while performing perturbative 
calculations. These problems come from the severe infra-red divergences which plague finite 
temperature Green's functions, in particular in gauge theories. To this aim a technique for 
the resummation of the so-called "hard thermal loops" was proposed and developed [ f2~8|j . 
However, this method is not successful for scales in which non-perturbative physics (such 
as the gluon magnetic mass in QCD) become relevant. The key problem is that ordinary 
perturbation theory does not clearly separate the various scales in the game, and the 
diagrammatic expansion has to be reorganized. 

We believe that the Wilson RG approach should be very helpful in this sense. 

In this paper we have considered the Wilson RG formulation of a RT thermal field 
theory. The main idea is to consider the thermal interaction between a quantum field and 
a thermal bath as an effective interaction, namely to regard the frequency modes above a 
certain scale A as effective interactions for the low energy modes below A. Therefore the 
thermal fluctuations above A are integrated out and put in a 'Wilsonian effective action" , 
which satisfies an "exact" (in principle) evolution equation in A. We solved this flow 
equation for the effective potential in various approximations, providing a non-perturbative 
resummation of Feynman graphs. The approach is physically quite transparent, rigorous, 
and gives better numerical results with respect to usual resummed perturbation theory. 
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The Wilson method had already been applied to the same self-interacting scalar model in 



the Imaginary Time Matsubara formalism of thermal field theories \W, [T7[]. Our results 
agree with these previous analyses and are better in some cases. However, the main aim of 
this paper was to set up the general formalism and discuss some approximation methods, 
by using the well-studied scalar theory. 

Concerning the extension of this RT approach to gauge theories, we believe it is very 
promising, since the issue of gauge invariance can be considered in a clean way by using the 
CTP formalism described in the Appendix A. Moreover, if one is interested in the static 
quantities (such as effective potentials and thermal masses), the imaginary parts of the self- 
energies can be neglected and only on-shell modes contribute to the flow, as shown by (ffl). 
In this case, by imposing the thermal boundary conditions only on the physical degrees of 



freedom, as suggested in |£9] , one finds that the thermal on-shell part of the propagators 
is gauge- invariant. Therefore the RG evolution equations obtained in this way provide a 
gauge-invariant resummation. The consequences of this observation will be considered in 
a separate paper. 

The RT formalism is absolutely necessary in order to deal with non-equilibrium phenom- 
ena. We hope this paper is a first step towards a consistent formulation of coarse-graining 
in non-equilibrium field theory. As far as the problem of thermalization is concerned, it 
should be possible to derive a Boltzmann equation, at least for a quasi-equilibrium or quasi- 
stationary system, and give a precise physical meaning to the IR cutoff A, which could be 
interpreted as the thermalization scale, A(£), to be determined dynamically. With such a 
systematic treatment, one should be able to compute corrections to the Boltzmann equa- 
tion, at least in some regimes. More generally, the dynamical elimination of the degrees of 
freedom which are irrelevant at a certain scale, would permit to derive an effective theory 
(for instance a kinetic theory |BD[) from the underlying fundamental field theory. These 



points are currently being explored. 
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A The Closed Time Path formalism 



In this Appendix, we give a short review of the CTP formalism [Q, and discuss its appli- 
cation to the formulation of the Wilson RG presented in this paper. We will follow the 
presentation of refs. || and || (see also the review [[31]). 

This technique was introduced to describe the evolution of an arbitrary initial state, 
described by a density matrix p. 
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In the usual path integral approach to scattering theory, we are interested to the tran- 
sition amplitudes between a \in) and a \out) vacuum state in the presence of an external 
source. We now face a different physical problem. We want to follow the out-of-equilibrium 
time evolution of an interacting field, without knowing its state in the far future. More- 
over, we look for a real and causal evolution equation. For a Hermitian field operator, the 
diagonal matrix element (in\in) is real and the off-diagonal ones (out\in) are complex. 

Therefore the basic idea of the CTP method is to start from a diagonal (in\in) matrix 
element at a given time t — and insert a complete set of states at a later time t'\ 

Z[J+,J-] = J[diP](in,0\7p,t')j_(iP,t'\in,Q)j + , (64) 

where J + and J_ are the two sources for the propagation forward and backward in time, 
respectively. Notice that one puts different sources for the two transition amplitudes, since 
one is mainly interested in following the forward time evolution. By introducing the density 
matrix p of the initial state and setting t' = +00 one gets 

Z[J+, J_,p] = J [d(p][d<f/][di/>](tp,0\T*exp (-i J J_$) |^,+oo) 
x(</>,+oo|Texp fi J J + $j \cp',0){(p',0\p\cp,0) , 

where T and T* are the time and anti-time ordering operators, respectively. We see the 
origin of the name CTP. This path-integral is performed on a contour in the complex time 
plane running from t = to t = 00 and then back to t = 0. One can then introduce 
a standard path integral representation for the two matrix elements corresponding to the 
transitions from the state at t — to the state at t = t' and viceversa: 

Z[J+, J_,p] = J [dcf> + }[d<f>„)(<f> + ,0\p\^,0)expi{S[<f> + ) + - S*[4>-] - J_0_} . (65) 

It is clear how the +, — fields and sources correspond to the ones labelled with 1, 2 in the 
paper, respectively. Therefore we introduce the vectors 

a =(0+,0_), J a = (J + ,J_). (66) 

As far as the initial state is concerned, one can always decompose the density matrix 
element at t = in the following way: 

(<h,0\p\<h,0)=expihc + J K^a + l J K ah <p a <p b + ...} . (67) 

Therefore the information on the initial state is all contained in the non-local sources K, 
which are obviously concentrated at t — 0. The various K's are just boundary conditions 
in time for the corresponding Green functions. In this way a perturbative expansion may 
be constructed, with the usual Feynman rules. 

We now specify the form of the initial state. For an initial state of a real free particle 
in thermal equilibrium, p ~ exp[— (3 J d 3 k Ukalak] and the matrix element is ||: 

(0 1 ,O|p|0 2 ,O)~exp-i J 51 sin ^ Ul(k) + 4>\{k)) cosh(/^ fc ) - 2Mk)Mk)] ■ (68) 
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We see that the effect of this initial state is just to add a correction to the free propagators 
of the theory. In this way the free propagator of the RT thermal field theory is obtained. 
For a interacting particle, the matrix element in (|67| ) cannot be computed exactly, but 
the Wick theorem may be used to perform the usual perturbative calculations of the RT 
equilibrium theory. In this way we find the same result as we would have found by applying 
the well-known property that, for a pure thermal state, the matrix element O|p|02, 0) 
admits a path integral representation over the path going from t — along the imaginary 
time axis to the time t = —i(3. It may also be shown that indeed there is a certain freedom 
in the choice of the path to go from t — to t — —i/3, so that all the various formulations 
of the equilibrium thermal field theory are obtained. 

The density matrix relevant to this paper, eq. @, is obtained by introducing a momentum- 
dependent temperature 





for 


|jfe| 


> A 




< 








(69) 


+oo 


for 


|jfc| 


< A 





As discussed in 0, if the density matrix has this form all the initial n-point correlations, 
K in (|5"7|), vanish for n > 2, and the result for the density matrix element ([IT]) is obtained 
from eq. fl68| ) after the substitution f3u k — > f3 k ^. 

As a consequence, the only modification with respect to the usual equilibrium RT 
formalism will consist in replacing the Bose-Einstein distribution function in the free prop- 
agators, with the cut off one defined in ([□]). 



B Cancellation of pinch singularities 

In this appendix we discuss the mechanism of cancellation of the pinch singularities which 
appear in the matrix kernel K\ ([36]), which we rewrite in the following way (we suppress 
the index A and restore the field indices 1,2 for sake of notation) 

1^ = 1 tri^rfj,. (70) 

We will be able to rewrite this equation in a way which is manifestly free of the pinch-like 
singularities A™Aq™ appearing in eq. 



First of all recall that the full matrix propagator G\ = [D]^ 1 + Sa] _1 is free of pinch 
singularities due to the matrix structure of D\ and Sa J?J, and so it is its A-derivative. 
Then the origin of the pinch singularities in the kernel of the RG flow equation ( [70]) is clear 
if one realizes that the definition ( |36| ) corresponds to the A-derivative of the full propagator 
taking into account only the A-dependence coming from the free propagator D\. This 
destroys the matrix structure of G\. In order to recover it (and therefore to cancel the 
pinch singularities), from the vertex in ( |T0"D there should come a contribution like 

-G A 
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which, added to the kernel (0), would give the total (pinch-free) A-derivative of G\. 

Indeed this happens in the following way. The A-derivative of the matrix self-energy is 
found by derivating with respect to <pi, <pj the evolution equation (j32|) 

Then one has the identity 

i^lij = iTxJj + f [^ ab G aa rGbb'[j^a'lrij ~ ^a'li^cd^b'dj] 5 (72) 

~ (3) 

where T k [■ is a vertex which is two -'particle-irreducible if one tries to separate the external 
fields ij from the external field k. By inserting this equation in ( |f0|) and using (|71|) one has 

9 r« - i tr 6 9 r \ 

which does not exhibit pinch singularities. We can summarize this mechanism of cancella- 
tion as follows: in the product K a i,T\^ h the contributions corresponding to the A-derivative 
of the matrix self-energy come from the parts of T U j in which the external fields ij are 
two-particle-reducible with respect to the external field 1 (the reader may convince himself 
of this by drawing the corresponding diagrams). These two-particle-reducible parts are 
the last term in the r.h.s. of eq. (|72|) . Then the A-derivative of the matrix self-energy 
cooperates with the kernel K\ to give the A-derivative of the full propagator, which is free 
of pinch singularities. 
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